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^ ■ We revisit unitary representation of centrally extended psu(2|2) excitation superalgebra. We 
find most generally that 'pseudo-momentum' , not lattice momentum, diagonalizes spin chain 
^ '. Hamiltonian and leads to generalized dynamic spin chain. All known results point to lattice mo- 
■ mentum diagonalization for =4 super Yang-Mills theory. Having different interacting struc- 
ture, we ask if fA/^ = 6 superconformal Chern-Simons theory provides an example of pseudo- 
momentum diagonalization. For S0(6) sector, we study maximal shuffling and next-to-maximal 
shuffling terms in the dilatation operator and compare them with results expected from psu(2|2) 
superalgebbra and integrability. At two loops, we rederive maximal shuffling term (3-site) and 
find perfect agreement with known results. At four loops, we first find absence of next-to- 
maximal shuffling term (4-site), in agreement with prediction based on integrability. We next 
extract maximal shuffling term (5-site), the most relevant term for checking the possibility of 
pseudo-momentum diagonalization. Curiously, we find that result agrees with integraility pre- 
diction based on lattice momentum, as in iA^^ = 4 super Yang-Mills theory. Consistency of our 
results is fully ensured by checks of renormalizability up to six loops. 



1 Introduction 



The holographic duality between Type IIA string theory on AdS4 x CP^ and (2+l)-dimensional 
9{, = 6 superconformal Chern-Simons theory, discovered by Aharony, Bergman, Jafferis and 
Maldacena (ABJM) yj, opened a new avenue for exploring the AdS/CFT correspondence |l2l. 
The (2+l)-dimensional superconformal field theory admits Lagrangian formulation, hereafter 
referred as ABJM theory, so a systematic study between bulk fields and boundary operators is 
feasible much to the same extent as the holographic duality between Type IIB string theory on 
AdSs X §^ and (3+l)-dimensional 9{, = 4 super Yang-Mills (SYM) theory. The ABJM theory 
depends on two parameters, level k of the Chern-Simons term and rank of the gauge group. 
The planar limit, — t- oo, of the ABJM theory corresponds to non-interacting limit, gs 0, of 
the Type IIA string theory. In this limit, the ABJM theory is organized in perturbation theory at 
weak 't Hooft coupling regime X={N/k)<til, while the Type IIA string theory is organized in 
perturbation theory at strong 't Hooft coupling regime y/% ^ 1 . 

A very interesting and important aspect of the AdS4-ABJM holographic duality is the 
prospect of 'integrability', the feature investigated thoroughly in the AdSs— SYM holographic 
duality [|3l- [[T6l . At strong 't Hooft coupling regime. Type IIA string worldsheet dynamics was 
found integrable at leading order in 1/ y/X ifTTl [TSl . At weak 't Hooft coupling regime, S0(6) 
sector of the ABJM dilatation operator was found integrable at 2 loops [fT9l[T8l . Extension to 
full OSp(6|4) dilatation operator and Bethe ansatz equations were proposed to all orders and 
numerous consistency checks of the proposal were studied [lIOll - ETl . It was found that the 
excitation spectrum in the AdS4-ABJM theory is organized by a direct sum of p5u(2|2) super- 
algebra. Interestingly, a direct product of the same superalgebra also featured in the excitation 
spectrum in the AdSs-SYM theory [|T5l . lfT6l . These results all point toward exact 'quantum 
integrability' of the ABJM theory, but a direct analytic proof would be highly desirable. 

The purpose of this paper is to study the dilatation operator of the ABJM theory beyond 
leading order. Our motivations are primarily twofold. First, we wish to understand precise dy- 
namic nature of the alternating spin chain in this theory. In the case of =4 super Yang-Mills 
theory, dilatation operator was described by a novel spin chain whose length changes dynami- 
cally [i7j. This novelty was in fact a simple consequence of psu(2, 2|4) superconformal algebra. 
For the osp(6|4) superconformal algebra of the ABJM theory, we also expect emergence of 
dynamic spin chain. Magnon excitations in both theories are known to be described by product 
or sum of centrally-extended psu(2|2). On the other hand, field contents and their interactions 
as well as supersymmetry preserved by the spin chain ground-state differ for the two theories. 
Therefore, precise dynamic nature of the spin chain could also differ each other. Second, we 
would like to test quantum integrablity of the operator contents in this theory. While there were 
several indications from both weak and strong coupling regime, there was no rigorous proof yet. 
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We address this problem by comparing spin chain Hamiltonian expected from integrabiUty and 
weak coupUng expansion of it with the dilatation operator computed directly from the ABJM 
theory at lower but beyond leading orders in perturbation theory. 

We begin in section 2 with centrally extended psu(2|2) superalgebra, the superalgebra that 
features both the ABJM and the fA^ = 4 SYM theories. Re-examining the previous study care- 
fully, we find that off-shell magnon excitations, described by unitary representations of the 
superalgebra, is in general labeled by a 'pseudo-momentum' which in general differs from the 
lattice momentum of the spin chain. Functional relation of the 'pseudo-momentum' to the lat- 
tice momentum depends on the theory under consideration, on the amount of supersymmetry 
and order of perturbative interactions. In turn, the 'pseudo-momentum' different from the lat- 
tice momentum lead to a generalized dynamical spin chain in that the exchange algebra of the 
magnon excitation is generalized order by order in perturbation theory. 

In section 3, focusing on a magnon in SU(2)®SU(2) sector, we compare predictions of 
quantum integrability to the spectrum of operator contents, equivalently, to general structure of 
the spin chain Hamiltonian at higher orders in perturbation theory. We focus in the dilatation 
operator on so-called 'maximal shuffling' terms - terms that exchange spins at furthest sites - 
and 'next-to-maximal shuffling' terms - terms that exchange spins one less than furthest sites - 
at each order, since these terms are the cleanest to compute yet highly nontrivial in structure. We 
explain diagrammatic origin of these terms up to four loops. At two loops, we rederive maximal 
shuffling term (3-sites) and find perfect agreement with known results. At four loops, we first 
find absence of next-to-maximal shuffling term (4-sites), in agreement with integrability. We 
then extract maximal shuffling term (5 -sites). 

In section 4, we explain formalism for computing anomalous dimensions from two-point 
function of single-trace operators in perturbation theory and hence the quantum dilatation oper- 
ator. We also explain general structure of the dilatation operators and how the maximal shuffling 
terms of length 2£+l and next-to-maximal shuffling terms of order 2£ arise diagrammatically at 
2£-th order in perturbation theory. In sections 5 and 6, we compute quantum dilatation operator 
in dimensional regularization. In section 5, we first compute two-loop contribution to two-point 
function and reproduce the known result of the dilatation operator. We then compute four-loop 
contribution. Combining both contributions, we extract next-to-maximal shuffling terms at four 
loops. If quantum integrability holds, result from section 3 indicates these terms should be 
absent. We indeed find that bosonic and fermionic loop contributions cancel each other - a 
clear indication that integrability holds. In section 6, we focus on maximal shuffling terms at 
four loops, where only bosonic diagrams contribute. Again, we find the result agrees with the 
prediction of integrability based on lattice momentum in section 3. Internal consistency of our 
results is fully ensured by checks of renormalizability up to six loops. 
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Along with = 4 SYM theory, our result escalates a puzzle why the dilatation oprators are 
always diagonalized by lattice momentum eigenstates, not by pseudo-momentum eigenstates. 
In section 7, we discuss implication of this to other issues pertinent to the ABJM theory. We 
also discuss various implications of our results. In Appendix A, we present psu(2|2) S-matrices 
in the pseudo-momentum basis. Appendix B is relegated to technical details of precision nu- 
merical evaluation of irreducible Feynman diagrams, whose results were used in sections 5 and 
6. 

2 Off-Shell ^eil(2 1 2) and Generalized Dynamic Spin Chain 

In the planar limit, we consider a single trace operator of the ABJM scalar fields. In ABJM 
theory, the scalar fields Y^y] [I = 1,2,3,4) transform as (N,N;4) and (N,N;4) under the 
U(A'^) X \^{N) gauge group and the SU(4) R- symmetry group. Consider in the planar limit 
single-trace operators . By gauge invariance, they form an alternating spin chain, labeled by 
SU(4) quantum numbers: 

[h hhhh--, hi] = Tr ( y^' yI yI y ^5 . . . y^t^ ) . (2. i ) 

In the limit L — t- 0°, choose the convention that Y' fields of 4 are in the odd sites of the spin chain 
and y/ of 4 in the even sites. By group theory, 4 at odd sites and 4 at even sites interact only 
through 5y. 

Rename {Y^ ,Y^ ,Y^ ,Y^) as (Ai,A2,52,5j) where Afl,^^ (a,b = 1,2) transform under the 
SU(2)^ and the SU(2)b subgroups of SU(4) ETTl . We focus on a restricted set of single-trace 
operators where Aa and field are placed at odd/ even sites, respectively: 

0[aia2 '336(4 6(5 ■■■] =Tr{AaiBa2Aa^Ba^Aa^ ■••). (2.2) 

A given operator of this type is mappable to a given spin chain state | ai a2 6(3 ■ • ■) . We shall take 
the ground state as 

|0) = |1111---) ^ 0[llll---]=Tr{AiBiAiBiAiBi ■■■). (2.3) 

For this type of spin chains, there is no interaction between 4 and 4 spins. Hence, the odd- 
site spin chain and the even- site spin chain behave independently. Then, there will be A-type 
magnons and 5-type magnons propagating independently without interactions between the two 
types 0. 

^In principle, interaction of double-exchange P„+i,n+3 Pn+2.n+4 could arise, where P,; denotes the permutation 
operator exchanging spins at the ;-th and the j-th sites. In this case, there will be interactions between these two 
types of magnons. 
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The magnon excitations are organized by off-shell psu(2|2) ©psu(2|2) superalgebra sym- 
metries acting on A- and 5-sites, respectively. The (2|2) excitation multiplets are (A2, ^2 l^i 5 ^2) 
and (52, 1^5,^1). Therefore, magnons of A- and 5-types behave identical. From now on, we 
shall focus on A-type magnons and denote (A2,52|^i,^2) as 4> = collectively. 
One expects a magnon excitation is an eigenstate of lattice momentum of the A-type sites in the 
alternating spin chain: 

L-l 

|p)=£e'"/'|...^...). (2.4) 

n=0 

Here, the ellipses denote spin configurations of the ground state, where A 1 or 5 1 in the remaining 
odd or even sites. Below, we shall suppress denoting these background spins explicit. 

The off-shell psu(2|2) superalgebra of the excitation symmetry is spanned by the two su(2) 
rotation generators 2^^, the supersymmetry generator 0^ and the superconformal gener- 
ator 6^. The off-shell configuration is characterized by 5/(2, M) central charges [[TSl . 
Their (anti)commutators are given by [fTSl 

{£1^, nj} = £«P£«,i?, {e^, eg} = e^^e^'Si* . (2.5) 

The central charges C is related to the energy by £ = C, while M* introduced at off-shell are 
related to the momentum. On a state of fundamental representation, the generators act as 

= dim - , ^W) = Spi^") - Iw^ ^^.6) 

and as 

6^1/) = JSSr). (2.7) 

The function G(A+) is not arbitrary. Note first that A+ inserts the ground-state spin Ai into the 
odd site of the chain while A ^ removes one ground-state spinAi. Acting on an eigenstate (|2.4I) 
of the lattice momentum p, we have 

G(A=*=)=A± where \A^^) = e^'P\^A^) . (2.8) 
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This is the defining relation of the new marker function G(A ). 

Closure of the superalgebra on fundamental representation leads to the 'shortening condi- 
tion' 



ad-bc=l. (2.9) 

One also finds the central charges of total momentum yield 

^ \ (^) ^ab\^G{A+)) 

= cd\^G{A-)) (2.10) 

and the central charge of total dilatation energy yields 

='^{ad + bc)\^). (2.11) 

Quantum mechanically, functional form of G(A^) in (|2.7I) are subject to radiative correc- 
tions. These corrections can be extracted straightforwardly, for example, from explicit deriva- 
tion of quantum dilatation operator and utilization of the relations (|2.10l 12.1 II) or from oper- 
ator product expansions between fermionic charges 0, & and a pair of single-trace operators 
carrying bosonic and fermionic excitations, respectively. After radiative corrections are taken 
into account, the full-fledged quantum marker function G{A^) ought to solve the 'shortening 
condition' (12.91 ). It is an interesting open problem to find a complete representation theoretic 
classification of G{A^) as an exact, nonperturbative function of X. 

Here, we content ourselves for a particular ansatz of G(A^) motivated by perturbation the- 
ory. Introduce a notion of 'pseudo-momentum' P by the exchange algebra: 

P:=P{p) where G{A^)^ := e^'^^P^^G{A^) . (2.12) 

Physically, we expect that P{p) is an odd function of p and periodic with periodicity In. This 
means that the pseudo-momentum P is expandable in a power series of the lattice momentum p 

as 

e'-f = G{e'P) = e'Pe'^n=i *2„(^) sinnp _ ^2. 13) 

From general structure of the perturbation theory that dictates hopping of an excitation spin 
over lattice sites, we expect that binQ^) is further expandable as 

oo 

^2„=£^2^,2a'^. (2.14) 
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To obey the on-shell conditions ^ = ^* = 0, one needs to excite two or more spins. Consider 
a spin chain with n multiple excitation spins that are asymptotically separated and carry 'pseudo- 
momentum' Pi, - ■ ■ ,Pn. The momentum central charges now read 

^*|4)i4)2 4>3...4,n^ =K*\^^^^^^- -^"G{A )) (2.15) 

with 

n n n n 

K=Y,akbk n K* = Y,Ckdk n (2-16) 

k=l m=k+l k=l m=k+l 

Both should vanish on any on-shell configuration that satisfies 

n 

Yle'^^l (2.17) 

In terms of the lattice momentum, this in general puts the spin chain to obey a version of twisted 
boundary condition if closed |§. Therefore, one might opt to consider the pseudo-momentum 
defined only in asymptotic limit. Our point is to stress that the condition at quantum level ought 
to be the one for the 'pseudo-momentum' P, not for the lattice momentum p. As shown in [flSl . 
a unique local solution is given by 

ukbk = a(?i)(e-'^* - 1) ; ckdk = |3(?i)(e'^* - 1) . (2.18) 

Note that a and |3 are independent of the lattice site k and the lattice momenta p^s - they are 
parameters common to all excitations. 

A convenient parametrization of a, b, c and d are in terms of elliptic variables and 

a = ./m Y, b = -^)J^{X+ -X-) 

c=v^^, d = -^h{}C}-{X^-X-). (2.19) 
With positive definite h{X), defining relations are 

^ = e'^; (X+)*=X-, (2.20) 

X 



while the shortening condition (|2.9I) reads: 



+ — = (2.21) 

X+ X h{X) ^ ' 



We thank N. Beisert for useful correspondences on issues related to this point. 
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Explicitly, 




Unitarity of the representation demands that, modulo a complex phase, y = 
and / = 1. In terms of these variables, the dilatation energy spectrum reads 




Classically, the 'pseudo-momentum' P is reduced to the lattice momentum p, and the elliptic 
variables are reduced to where = e'^. 

As in ifTSl , S-matrices between a pair of magnon excitations is determinable by the psu(2|2) 
superalgebra. For completeness, we tabulate them in appendix A. Following [1281 , overall scalar 
phase-factor is also determinable by imposing crossing relations. These S-matrices satisfy uni- 
tarity, Yang-Baxter equations and crossing symmetries. From the S-matrices, one can also 
construct full-fledged Bethe ansatz equations (BAE). We stress that all these conclusions are 
most transparent when the 'pseudo-momentum' P, equivalently, the elliptic variables X^ are 
used instead of the lattice momentum p, equivalently, . 

A remark is in order. Deformation of integrable spin chain was considered in f29\, in which 
the coupling parameter X is replaced by a set of four deformation parameters. In particular, one 
of these parameters acts to deform {x± + /x±) — > jc± + L;i=3 oCn/ {^±)" and hence the energy 
spectrum as well. It was shown that these deformations are possible while retaining global 
symmetries of the spin chain intact. Such a deformation appears closely related to the map 
from the lattice momentum to the pseudo-momentum introduced above. It would be extremely 
interesting to understand possible relation better. 

3 Spin Chain Hamiltonian from Integrability 

The magnon spectrum expected from the psu(2|2) excitation symmetry and the integrability is 
of the form 



It was shown that the energy spectrum receives perturbative corrections only from even-loop 
orders. So, the interpolating function Ji^iX) is parametrizable as 




(3.1) 



with 



(3.2) 
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This is the spectrum for the single magnon for each SU(2) sector. Expanding the energy spec- 
trum at weak coupling as a function of lattice momentum p, 

E = l + 4/z2(X)sin2|-16/z4(X)sin4| + ... 

= Q) + (4sin2|)x2+(-16sin4|+4/.2sin2|)xV-- . (3.3) 

As is well-known from planar perturbation theory, 2£-th order contribution to the energy gives 
rise to lattice shuffling up to 2£ consecutive sites. At 2-loop order, from e^'^, the maximal 
lattice shuffling is for 2 sites. At 4-loop order, from e^'^'P, the maximal lattice shuffling is for 4 
sites, etc. We shall refer to these as 'maximal shuffling' interactions. 

At this stage, as shown in the previous section, one needs to bear in mind of the possi- 
bility that dynamic spin chain is diagonalized in the pseudo-momentum P, not in the lattice 
momentum p. In this case, we need to replace the lattice momentum p in (13.11) by the pseudo- 
momentum P. As the pseudo-momentum is defined as a function of the lattice momentum in 
perturbation theory, P = P{p), the energy spectrum is still expandable as double series of the 't 
Hooft coupling and the lattice hopping. Schematically, it takes the form 

i=0 n=0 ^ 




(3.4) 



where eo.o = \ counts the classical scaling dimension and ^2^.0 = (£ > 1) is fixed by the 
supersymmetry of the ground-state. The coefficient ^2.2 = 4 was computed previously from 
explicit computation of quantum dilatation operator at two loop order lfT9llT8l . 

The coefficients e2£,2n provide a complete information concerning spectra of the quantum 
dilatation operator, E{p,X). More precisely, from these coefficients, one can determine both 
the generalized marker function G{e'P) and the interpolating coupling function h^Ck) as the two 
functions are mutually independent. Our assertion is based on the observation that, at each order 
in perturbation theory, the coefficients e2^,2n in the spectra are in one-to-one correspondence 
with the two sets of coefficients (/i2f,^2^,2n)- 

The corresponding integrable Hamiltonian at each order is well-known. Consider the odd- 
site SU(2) chain. The tree-level part of the Hamiltonian counting the classical scaling dimension 
is 

1 

Ho = -Y,^- (3.5) 

n=0 
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The 2-loop part of the Hamiltonian is simply given by [fT9l [T8l 

^2=i;(i-p«,„,,.,„,3): (3-6) 

11=0 

where 2L is the number of the total sites and P denotes exchange operator. We consider the 
asymptotic spin chain where L goes to infinity. The 4-loop part of the Hamiltonian can be 
identified as 

^4=^1(1- P«2„+,.2„+3 ) + I [ (P«2„+l«2„+5 - 1 )+ 4(1 - P«,„,,«,„,3 )] , (3.7) 

where we have used the result of the 9{, = 4 super Yang-Mills theory [31 171 HI. In order to 
demonstrate that the energy spectrum (13.41 ) follows from these Hamiltonian, we use the mo- 
mentum eigenstate of a single A-magnon excited on the odd-site SU(2) chain: 

\p)A=Y^e"^P\...^...) (3.8) 

n=0 

in the background of Ais on the ellipses. Likewise, for the even-site SU(2) chain, the momentum 
eigenstate of a single 5-magnon is 

\p)b= £e""^'|.--"Br...) (3.9) 

m=l 

in the background of Bis on the ellipses. Therefore, there are four kinds of excitation states: 

|0)a®|0)b, |p)a®|0)b, |0)a®|/)b, \p)a^\p')b. (3.10) 

The second and the fourth states are the momentum p states of the odd-site SU(2) chain. Below, 
we shall focus on the odd-site chain because the odd-site SU(2) chain and even-site SU(2) chain 
behave independently for the above quartet of states. 
One easily checks that 

2 P I \ 

^^2|;?)a = 4sm - \p}a 

Ha\p)a = (^64,2 sin^ ^ ^4^4 sin"^ (3.11) 

etc. In the Hamiltonian (13.71) . the maximal shuffling terms e2i,2i are of particular interest. For 
example, at 4-loop order, (13.31 ) indicates that 64 4 equals to —16. In perturbation theory, 64 4 in 
H/^ can be extracted by evaluating anomalous dimensions of a single A-magnon up to 4-loop 
orders. Since the value ^4,4 = —16 is a consequence of psu(2|2) excitation symmetry and the 
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integrability, explicit computation would amount to a test of quantum integrability of the AB JM 
theory. 

In the following sections, we shall compute the maximal shuffling term ^4 4 and the next- 
to-maximal shuffling term 643 explicitly and compare with prediction of the psu(2|2) exci- 
tation symmetry and the integrability. For the coefficient 64^4 of the maximal shuffling term, 
one needs to compute the diagram in Fig. [T] and its conjugate. Diagrammatically, there can 




Figure 1 : Maximal shuffling: 4-loop interactions involving five sites. 



also arise next-to-maximal shuffling terms. To compute their coefficient ^4,3, one needs to 
compute the diagrams in Fig. |2l and Fig. [3l They give rise to dilatation operators of the form 
^„ (I — Pn+i ,„+3 ) (I — P/i+2,n+4 ) • In case one type- A and one type-B magnons are excited simul- 
taneously, this operators represent nontrivial interactions between them. This operator, how- 
ever, is not diagonalized by the state | p)a ® | p')b- Absence of such interactions amounts to the 
statement that 643 =0. 




Figure 2: Next-to-maximal shuffling: 4-loop interaction diagram involving four sites. 



Figure 3: Next-to-maximal shuffling: 4-loop interaction diagram involving four sites and fermion loop. 
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4 Quantum Dilatation Operator 



O{x)::O{0):), = ——^^-—— (4.1) 



4.1 Formalism 

We use the anomalous dimension matrix that is defined by the two-point correlation function: 

(;c2)^-2^(;c2)Y(e) 

where we use the dimensional regularization to control the ultraviolet divergences. This is 
related to the collections A„ of Feynman diagrams at each order defined by 

(: 0{x) : : 0(0) :)e = (4)2ie-Yein(x2A2(e)) 

:=(/e)2^exp[ln(l+A2X2+A4X4^---)] . (4.2) 

Here, 4 denotes the dimensionally regularized Euclidean scalar propagator in the position 
space: 



/ = / ^ 1 e-P- = ^^^-^^ -J— (4 3) 



with 2(0 = 3 - £. 

From the above, we extract lower loop contributions to the dilatation operator in the AB JM 
theory § by 

//2 = -lime^oeA2 
H4 = -lime^o 2£ (^A4 - ^aI 

He = -lime^o 3£ (as - ^ (A2A4 + A4A2) + ^A^^ . (4.4) 

Here, the overall factor i in H2£ arises from extracting the dilatation operator as the charge 
generating variation of lnx2 in the ABJM theory. The leading singularity in A2 starts from 
the order 0(£^), so the leading term in the Laurent expansion of A2 contributes to H2. For 
the case of A4, the leading singularity in general starts at the 0(£^2) order. The coefficient of 
this leading singularity in A4 — ^A^ ought to vanish for renormalizability of the theory under 
consideration. Therefore, the leading singularity starts again from 0{e^^) terms. For the case 
of Ag, the leading singularity starts at the 0(£~^) order. Again, the coefficient of this leading 
singularity of 0{e^^) and of the next-to-leading singularity 0(£^2) vanish for renormalizability. 
This pattern continues to all higher order contributions to the dilatation operator H. 



^For perturbative computations of the dilation operator in 9{, =4- SYM theory in close parallel to the Feynman 



diagrammatics we take, see ll30l . 
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4.2 Maximal Shuffling and Next-To-Maximal Shuffling Interactions 

We now focus on the ABJM theory. The operator structure involved with Fig. [T]can be con- 
structed from the operator structure of the two-loop scalar sextet contribution: 

Oj23 = 21 - K„+i^„+2 - IK„+2,n+3 + ^^n+i,n+3'^n+l,n+2 + 2P,;+l,„+3lK„+2,n+3 " 4P„+l,„+3 (4.5) 

where (K)fj^ = 5/y5^^ refers to contraction operator. The matrix structure of Fig.[T]is then 

^^3450r23 = 41 - 8P„+l,„+3 - 8P„+3,n+5 + 16P„+3,n+5P«+l,n+3 + " " " (4.6) 

where the ellipses denotes the omission of the terms involving any types of contraction between 
the odd and the even sites. Similarly, there is also a diagram corresponding to the conjugate of 
Fig.m For this, the operator structure takes the form, 

OmOsAS = 41 - 8P„+i,„+3 - 8P„+3,n+5 + 16P„+L„+3P„+3,n+5 + " " " (4.7) 

Therefore the full operator structure becomes 

*^12345 — 123 '^345 +'^345*^123 

= 81 - 16P„+i,„+3 - 16P„+3,n+5 + 16P„+L„+3P.+3,n+5 + 16P.+3,„+5P„+l,.+3 + " " <4.8) 

If we put only Aa I Ba magnons to the odd- / even-sites of the chain, the above can be rewritten 

as 



0^2345 = 16P„+l,„+5-8I + -- - , 

where we have used the identity 

tla,a,a,tlb,b,b, = I + 5,^J5,^^5^^ + S^JS^^S^^ - P13 - P35 - P15 = . 



(4.9) 



(4.10) 



The same operator structure arises from —\a\, the second-order effects of the two-loop con- 
tribution to the dilatation operator. Hence, adding these two contributions, one obtains the 
coefficient ^4 4. 



Figure 4: Remaining 4 loop interaction diagram involving four sites 
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The operator structure of the diagrams in Fig. [2] involves 4-site interactions and can be 
extracted as 

^1234 — ^^123*^234 + *^234*^123 

= 32(P„+i,„+3 - I) (Pn+2,„+4 - I) + 8(2P„+i,„+3 + 2P„+3,,+5 - 31) + ■ ■ ■ .(4.1 1) 

Again, the same operator structure arises from —\a\ contribution to the dilatation operator. 
Fermions also contribute to the 4-site operators Pn+i,n+3lPn+2,;i+4- This involves fermion loop 
as depicted in Fig.[3]and the corresponding matrix structure is given by 

<5r234 = Of234 + 32I + ---. (4.12) 

Operators of the form (Pn+i,n+3 — I) (Pn+2,n+4 — I) contributes to ^4 3 in 7/4. As said, the 
p5u(2|2) excitation symmetry and the integrability predicts that ^4,3 = 0. In the following sec- 
tions, we shall confirm this explicitly — the contributions of Fig. |2l Fig. |3] and second-order 
effects of 2-loop contributions all cancel one another. 

There are many other diagrams that arise at 4-loops and contribute to lower shuffling op- 
erators. Though we shall not evaluate any of these terms in this paper, for completeness of 
our discussions, we shall list a class of relevant diagrams. The operator structure of Fig. |4]is 
proportional to 

er234 = 4P„+l,„+3 + 4P„+2,„+4 - 41 + ■ ■ • (4. 13) 

Therefore this 4-loop diagram is relevant to the computation of the coefficient £4,2- 



Figure 5: 4-loop diagram involving 5 sites and only K operators 





Figure 6: 4-loop diagram involving 5 sites and only K operators 
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Figure 7: 4-loop diagram involving 4 sites and only K operators 



>-< 



Figure 8: 4-loop diagram involving 4 sites and fermion loop 




Figure 9: 4-loop interaction diagram involving 3 sites 




Figure 10: Another 4-loop interaction diagram involving 3 sites 



There are many 4-loop diagrams involving only K operators. Such diagrams involving 5- 
sites are illustrated in Fig. [5] and Fig. [6] . Further diagrams involving 4-sites are illustrated in 
Fig. n\ and Fig. [8l The 4-loop diagrams contributing operators of 3-sites are proliferated. In 
Fig. m and Fig. \T0[ we illustrate two representative diagrams arising from sextet interaction 
vertices only. They contribute to the coefficient 64^2- 



5 3-Site and 4-Site Interactions 

In this section, we shall first compute contributions to 3-site operators from 2-loop diagrams 
and 4-site operators from 4-loop diagrams. As discussed in the last section, the second-order 
effect of the 2-loop diagrams will also contribute to H4 along with 4-loop diagrams. 
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5.1 2 loops and 3-site interactions 



For the computation of the 4-loop diagrams, we begin with the computation of the 2-loop and 
the 1-loop diagrams with a momentum flow P. We first compute the 2-loop contribution of 
Fig.dUusing the dimensional regularization. In this diagram, we suppress symmetry factor first 
and then reinstate it in the end. We then have 

d2> d% 1 1 1 (r((0-l))3r(3-2(0) 



(27r)2» {Inf^ (^p + g + p)2 (4;^)2a)r(3(o_3)(p2)3-2a) 

To get this result, we have used 

r d2> {4Ti)'^{P^)"+b-^ ^r{a + b- a))r((0 - a)r(a) - b) 
^ ~J (27r)2» (p2)a((p + p)2)^, - r{a)r{b)r{2(0-a-b) ' 

Similarly, the 1-loop contribution in Fig.[TT]is evaluated as 

1 (r(a)-l))2r(2-a)) 



LiiP) 



d^'^p 1 



(27r)2»;?2 (p + p)2 (47r)"r(2(o-2)(/'2)2-co- 



(5.1) 



(5.2) 



(5.3) 





Figure 11: Basic 2-loop and 1-loop diagrams with a momentum How 



Figure 12: 3-loop sextet interaction diagram involving 3 sites 

Based on these building blocks, we next compute the diagram of Fig. [I2l It is proportional 
to the operator ^^123- 

K2 = {hfA2 = (/e)3(x27l)3-2co ^2 £0^23 , (5-4) 

where a2 is given by 



-3/^2^x2(0-3 



a2 = {-{2nr)-ih)-'ix'n 



2 JP x 



(27r)2t« 

(27t)2(r(a)-l))3(r(3-2co))2r(5a)-6) 
(47r)3(r(3co-3))2r(6-4a)) 



4£ 



[l+£(l-l|/(l/2))+0(£2)]. 



(5.5) 
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Here, the polygamma function \|/(l/2) takes the value: 



¥(1/2) 



-C-21n2: 



1.963510026. 



(5.6) 



where C denotes the Euler's constant. For the Fourier transform, we used the formula, 



(f'^p 1 



r((o-a) 1 



(5.7) 



(271)2" (p2)a 4«7r"r(a) (x2)'^-« ' 

We see that the part of the 2-loop Hamiltonian H2 proportional to P„+i.;3+3 is solely coming 
from this sextet interaction contribution A2. Therefore, this part of H2 becomes 



H 



3s 



-lime 



1 



1 



^e^0ea2lCf23 = -£Or23=I 

This reproduces the previous result obtained in Refs. [|T9l[T8l 



(5.8) 



5.2 4-loops and absence of 4-site interactions 

We now turn to the coefficient 64 3 of the operator 0(234- The contributions come from the 4- 
loop diagram of Fig.[2]and also from —\a\, the second-order effect of the 2-loop contribution. 
If the integrability were to hold, contributions to ^4 3 ought to be zero. 

We first compute the 4-loop contribution of Fig. [21 The full contribution including the (/e)"^ 
part takes the form 



(5.9) 



To proceed further, we also need the following double integral: 



72 = /sf 1 , 3 - 2(0, 1 , 3 - 2(0, 2 - (o) 



(5.10) 



where 



/5(wi,W2,W3,W4,W5) 



^2co^ ^2co^ 



(471)2^^ (p2^~2co+li^lW, 



{2n)2^ (27r)2» (/t2)>^i(Z2)W2((/_p)2)w3((^_p)2)w4((^_/)2)w5 

The coefficient Z74 can then be computed as 



(.5.11) 



Z74 = 72 



(r(a)-i))^ 
44712 



r(3-2(o) 



r(2-(o)r(8(o-io) 



(5.12) 



_r(3(o-3)J r(2(o-2)r(io-7(o) ■ 

In the Laurent expansion in powers of £, the integral 72 starts with the 0(e*^) order. Therefore, 
it can be expanded as 



72 = |3i l+p2e + 0(£2 



(5.13) 
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Then Z74 has the expansion 



b4 



Pi 



128TO2 

where the polygamma function is given by 



l + (|32-5-v(l)-2v(l/2))£ + 0(£2 



-0.577215. 



(5.14) 



(5.15) 



On the other hand, the corresponding coefficient of 0^234 from — is again given by —jo^- 
The leading singularity of order 0{e^^) contribution from these two ought to cancel each other 
since the ABJM theory is renormalizable. We shall shortly show that the 4-site contribution of 
Fig. [3] starts with the order 0(e^^). This then the coefficient Pi to 



|3i = -47r. 



(5.16) 



We checked this numerically with high precision in the Appendix B. The corresponding dilata- 
tion operator is then given by 



H 



We found that 



-lime^o 2£ 



b4 - ^aj 



1 



EOf234 = -- |32-7-V(l) IOf234. (5-17) 



^2 = S+W) 



(5.18) 



by carrying out the integral numerically. Again, see Appendix B. Therefore, we find scalar loop 
contributions sum up to yield 



1 



1234 • 



(5.19) 



Let us now turn to the contribution from the fermion loop contribution. Fig. [3l There are no 
more terms that are proportional to 0(234- Since this diagram is of 0(e^^), we may extract its 
contribution to the dilatation operator directly by computing the diagram in Fig. [13] with zero 
momentum flow to the amputated external lines. The corresponding Feynman integral is given 



1 X 

_A 



Figure 13: 4-loop diagram with amputated external lines 
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by 

where dji? abbreviates for d^^pj [liCf-^, etc. Then, reinstating the symmetry factor, the contri- 
bution yields 

=F4-^- (271)4. l£Of234 (5.21) 

To compute F\, we perform /-integral first, p and q integral next, and finally ^-integral using the 
definition /s and G(fl, h). One finds 

^4=^^G(l,l)G(5-3to,l)[2G(2-to,l)G(4-2to,l)-75(l,l,l,l,l-to)]. (5.22) 



One can evaluate 



G(l,l)G(5-3o),l) = ^(l+0(£)) 

2G(2-co,l)G(4-2to,l) = — [l + 8(3 + \|/(l)) + 0(82)] . (5.23) 



In the Appendix B, we further find numerically that 



Hence, 



73 =/5(l, 1, 1, 1, 1 -00) = Y [1 +>l'(l)) + )] (5-24) 



^^ = S-^L'5r234 (5.25) 

and the fermion contribution the Hamiltonian Ha, reads 

^4 ^ = lime^o2£A^ = 1 £ Of234 • (5-26) 

Adding the bosonic and the fermionic contributions, the total 4-site contribution becomes 

Hf = hI + h{ = Y,{1\ + ---) , {5.21) 

and e4,3 = identically. We see that this conclusion fits exactly with the psu(2|2) excitation 
symmetry and the integrability. 
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6 5-Site Interactions 

We now turn to the 5-site interaction generated by the diagram in Fig. [T]0 This computes the 
coefficient 64^4. Below, the double integral (15.111) for some particular choice of (wi , ■ ■ ■ , W5) will 
play a central role. We need to evaluate the integral 



7i =75(2 -(0,3 -2(0, 2 -co, 3 -2(0, 1 ). 
The whole contribution of Fig. [H then takes the form 



(6.1) 



(6.2) 



where 



CI4 = Ji ■ 



(r(a)-i)) 
44712 



5 r 



r(2-(o)r(3-2(o) 
r(2(0-2)r(3(0-3) 



r(9(o-ii) 
r(i 1-8(0) 



(6.3) 



In the Laurent expansion in £, 7i starts with 0{e ^) . Therefore, the integral 7i is expandable 



as 



CXl 

Ji = — 

£ 



1 +a2e + 0(£' 



Then a4 is expanded as 



a4 



1 

On the other hand, has the expansion: 



3ai7r 

'64£2 



1 + (a2 - 10 - 3\k( l/2))£ + 0(£2) 



(6.4) 



(6.5) 



1 2 1 
2''2 32£2 



l + (2-2\|/(l/2))£ + C>(£^ 



(6.6) 



Again, the coefficients of 0{e ^) in a4 — must vanish by the renormalizability of the ABJM 
theory. This determines tti as 



371 



(6.7) 



By numerical integration with high precision, we found that tti agrees with — 2/(37r) with high 
accuracy. Again, see Appendix B. The corresponding contribution to H4 is then given by 

1 



Hi' = -lime^o 2£ 



a4 — —6(9 
, 2 . 



n 

12345 



a2-12-\|/(l/2) ^P„+i,„+5 + 



(6.8) 



^We are grateful to Joe Minahan pointing out a potential error to the final conclusion in this section concerning 
(in)compatibility with lattice mommentum. 
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The coefficient a2 is computed in the Appendix B as 



a2 = 13 + v(l/2). 



(6.9) 



The H^^ is thus given by 



— — ^lP«+i,;j+5 H • 



(6.10) 



Comparing this with the 5-site maximal shuffling term in H4, we deduce that 



64,4 = -16. 



(6.11) 



This agrees precisely with the prediction, 64,4 = —16. This prediction was based on psu(2|2) 
excitation symmetry and the integrability. On the other hand, the prediction was also based on a 
tacit assumption that the magnon excitation is diagonalized by lattice momentum eigenstate. Is 
it possible that the magnon excitation is actually diagonalized by pseudo-momentum eigenstate? 
We answer this question in next section. 

We close this section with a remark on internal consistency check of the above computations. 
It is also possible to extract the 6-loop contribution to the Hamiltonian, H(, in (14.41) . There, we 
first need to make sure that terms of order 0{e^^) and 0(e^^) in [A(, — ^(A4A2 +A2A4) + ^A^] 
vanish identically by the renormalizability of the ABJM theory. Cancelation of these terms 
provide a set of stringent consistency test for the computation of 04 since combinations of 04 
and a2 (as well as cubic of a2) ought to cancel against a^. The check is straightforward. We 
confirmed that these singular terms in £-expansion indeed cancel, ensuring that we computed 
a4 correctly. Details of this computation, along with explicit evaluation of H(„ will be reported 
in a separate paper. 

7 Interpretation and Discussions 

7.1 Interpretation 

In the last section, explicit computations showed that the 4-loop dilatation operator shows the 
structure of anticipated maximal and next-to-maximal shuffling term. Moreover, the coefficient 
of them matched with the value (13.11) dictated by the psu(2|2) excitation symmetry and inte- 
grability. In this subsection, we compare the result with consideration of pseudo-momentum 
discussed in section 2. We explain that, as in the 9l = 4 super Yang-Mills theory, the dilatation 
operator of the ABJM theory is diagonalized by the lattice momentum eigenstate, not by the 
pseudo-momentum eigenstate. 
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To show this, we revisit the energy spectrum (13.11 ) of a single SU(2) magnon but now as- 
suming that it is an eigenstate of in terms of the pseudo-momentum. Using the relation (12.131) 



oo oo 



P{p) = + L E ^2£,2« ^'^^ sinnp 



and the interpolating coupling function (13.21) 



(7.1) 



h^{X) = X^Y,h2e^ with 



hn := 1 



the energy spectrum 



l+4/.2(^)sin2M 



is now expanded in terms of the lattice momentum as 



1 



4sin2^ 



; 16 + 8^72,2) sin^ I + (4/12 + 8Z72,2) sin^ | 



This falls into the pattern we conjectured for the most general spectrum (|3.4I) : 



£X2^£.2.,2nSin2"f 

"-0 n=0 ^ 

^ j + (.2,2 sin^ + (e4,2 sin^ ^ + ^4,4 sin"^ ^'^^ ^' 



(7.2) 



(7.3) 



(7.4) 



(7.5) 



with 



.2,2= 4, .4,2=4/12 + 8/^2,2, .4,4 = -16 -8/72,2,- ••• 

One notes that shift from lattice momentum to pseudo-momentum affects the coefficient of 
the maximal shuffling term. From the last relation one learns that our computation (16.1 II) 
amounts /'2,2 = 0, viz. lattice momentum eigenstate diagonalizes the dilatation operator. 

It is easy to see that the identification works uniquely at each order in perturbation theory. 
The coefficients e2£,2n of conjectured terms in the energy spectrum (13.41) is in one-to-one cor- 
respondence with the two sets of coefficients {h2i,b2£,2n)- Since the latter two originate from 
totally different functions, the interpolating coupling function h^Ck) and the generalized marker 
function G(z), we see that the spectrum X) determines these two functions uniquely. 

^The first relation corresponds to the 2-loop results lfT9]| ifTSl . The second relation carries useful information 
regarding the interpolating function h{X) to the first nontrivial order in the weak coupling perturbation theory. 
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7.2 Discussions 



In this paper, we studied realization of psu(2|2) exciation symmetry and potential integrability 
in the ABJM theory. The theory is very different from the fA^ = 4 super Yang-Mills theory, 
only sharing these two aspects. Given this, we raised the possibility that eigenstates (not just 
eigenvlaues) diagonalizing the dilatation operator are not the lattice momentum basis but the 
more general pseudo-momentum basis whose functional form P{p) is subject to perturbative 
corrections. Our result shows the contrary: the ABJM theory is essentially the same the 9{, = 4 
SYM theory, at least in these aspects. 

Our result confirms that the ABJM theory is quantum integrable up to 4-loop order in weak 
coupling perturbation theory. We focused on the SU(2)xSU(2) inside S0(6) scalar sector and 
on maximal shuffling operators in the Hamiltonian, but we expect the integrability extends 
to the full OSp(6|4) and to all operators straightforwardly, as in [|20ll . We believe our result 
adds further evidence for exact quantum integrability of the ABJM theory. Below, we discuss 
implications of our results and issues that deserve further investigation. 

• Our consideration in section 2 indicates that notion of pseudo-momentum basis is ubiq- 
uitous for the off-shell psu(2|2) superalgebra [15]. This superalgebra features both the 
iV; = 4 SYM theory and the ABJM theory, but since contents and details of the two the- 
ories are very different, precise form of the pseudo-momentum P{p) would differ for the 
two theories. Nevertheless, result of this paper indicates that the dilatation operators of 
both theories are diagonalized by lattice momentum eigenbasis. The fact that change to 
pseudo-momentum basis did not actually take place in both theories poses a puzzle and 
additional hidden structure yet to be understood better. A possible checkpoint would be 
a direct computation of the off-shell central charges in weak coupling perturba- 
tion theory. By comparing the result with the local solution (|2.18|) . one may be able to 
understand why all ^2^,2n vanish in these two theories. 

• Our result also bears an implication to the interpolating function h^{X). The leading 
perturbative correction h2 in (|7.2I) is extractable, for instance, from 3-site shuffling term 
€4^2 in the 4-loop spectrum §. From (|7.6I) . one sees that 64,2 is directly related to if ^2,2 
vanishes and the dilatation operator is diagonalized by the lattice momentum eigenstates. 

• The notion of pseudo-momentum also calls for revisiting the interpretation of the gi- 
ant magnon at strong 't Hooft coupling limit for both AdSs ll32l and AdS4 [|27l cases.. 
In the light-cone gauge of the Lorentzian string worldsheet, it is always possible to 

^The computation is currently under progress by J. Minahan, O. Ohlson Sax and C. Sieg (private communica- 
tions). See also lISTl . 
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reparametrize each chiral worldsheet coordinates separately. This implies that there can 
in general be an arbitrariness in identifying the worldsheet momentum conjugate to the 
chiral worldsheet coordinate with the angular separation of the string configuration in the 
giant magnon. 

• Extension of our result to the parity-violating 9^ = 6 superconformal Chem-Simons the- 
ory with gauge group U(M)xU(A^) is straightforward. Two-loop Hamiltonian was con- 
structed in [|33l and was found integrable. Our result in this paper indicates that the 
integrability persists up to four loops: one only needs to replace to the geometric mean 
ofl?{M) and I'^iN). 

• By introducing a variety of D-branes wrapping supersymmetric cycles in CP^, one can 
construct open spin chains. For AdSs / SYM4 correspondence, several open spin chains 
including those ending on giant gravitons were found quantum integrable [|34l . The 
psu(2|2) superalgebra or appropriate subalgebra continued to act as the excitation sym- 
metry. Similarly, we expect that integrable open spin chains exist in ABJ(M) theories by 
adding D-branes and flavors [|35l . With concrete realization and identification of excita- 
tion symmetries therein, we expect that one can learn more not only about the quantum 
integrability but also aspects of the pseudo-momentum in the presence of boundaries. 

• Extension to higher loop orders would be highly desirable. It will teach us not only 
the magnon energy spectrum and quantum integrability but also other pertinent issues 
including wrapping interactions, structural (dis)similarity of the dilatation operator with 
the fA(; = 4 SYM theory etc. 

We are currently investigating these issues and intend to report further results in separate 
publications. 

Acknowledgement 

We are grateful to Niklas Beisert, David J. Gross, Joe Minahan, Didina Serban and Matthias 
Staudacher for very helpful discussions and correspondences. SJR thanks warm hospitality 
of the Kavli Institute for Theoretical Physics during this work. This work was supported in 
part by by Grants from the Korea National Science Foundation RO 1-2008-000- 10656-0, SRC- 
CQUeST-Rl 1-2005-021, 2008-313-C00175, 2005-084-C00003, 2009-008-0372, EU-FP Marie 
Curie Research & Training Networks HPRN-CT-2006-035863 (2009-06318) and the U.S. Na- 
tional Science Foundation under Grant No. PHY05-51 164 at KITP. 



23 



A p5u(2|2) S-matrices with Pseudo-Momentum 

In this section we record the SU (2|2) invariant S matrix. Only change is to replace of S 
matrix in [fTSl by which have more complicated momentum dependence. The S matrix 
readJzl 



where 



5l2|V?'(^^) =iS:i2|V2>l)+^12|(^2¥?') (A.l) 



'^2' — 



S _ ^0 ^2 -^f f 1 3 1 - 1/(^2"V) X^-X^ 

'^X^ ~X+\ 1 - l/iX^X^) X+-X, 



C12 



^0 2yiY2/(X-Xf) 
''^\-l/{X^X^) X^~X+ 



D12 = -5i2 

^12 = -^?2(l-2 



1-1/(X2+Xf) x^-x^ 



1-1/{X+X+) x^-x^ 



12 



d)T2 xr-xr 

A12 = 017 X 

^12 = ^?2 ^'-I^+ • (A.2) 

The phase 5]*2 satisfying 5'^25'2i = 1 will be specified below. The above is for the odd chain 
S matrices. For the even chain one can have 2|2 excitations transforming under the 5'm(2|2) 



^More precisely one needs to include the correction in Ref. lfT4ll . 
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symmetriy. The S matrix for the N = 6 Chern-Simons theory is consisting of 

5^^ = 5^" = 5?2S. (A.3) 

where S denotes the matrix part of the above 5 matrix without the phase • One also finds the 
crossing symmetric phases by 

n l-l/(X+X-) , 
12 1-1/(X2-X+) ^ 

^~?2 = ^+I^|- <^(^b^2) (A.4) 



where QiXx^X^) is the dressing phase found in [fT3l . 

The full 5 matrices above satisfy the unitarity condition 



and the Yang-Baxter equation 



s{Ax=l (A.5) 



We note that these structures relies only on local exchange algebras, independent global bound- 
ary conditions such as (|2.17l) 



B Numerical Evaluation of the integrals /i, J 2 and 

In this appendix, we evaluate the two-loop integrals /5(wl,w2,w3,w4,w5) in (15.1 II) . It is 
difficult to find an analytically closed form of them. However there exist several packages 
with which we can numerically evaluate them. It turns out that the Mathematica packages, 
MB[37l[38| and AMBRE[[39l are useful. Since these packages utilize a method based on the 
Mellin-Bames representation of integrals, we cast the integrals into the form: 

' dz\ f (iz2 r(-zi)r(a)-w25 -zi)r(a)- wi 

r(wi -zi] 
r(-z2)r(a)- W35 -z2)r((o- W4+Z2) 



/5(wl,w2,w3,w4,w5) = / — / — — ^ (B.l) 

J 2ni J 2ni I wi — zi 



X 



r(w4-Z2) 

r(-(0 + Wi4 -Zl -Z2)r(-(0 + W235 +21 +Z2)'r{ws+Zl + Z2) 



r(w2)r(w3)r(w5)r(w235)r(2(0 - W235)r(2a) - W^+Zl +Z2) ' 

where W235 denotes W2 + W3+W5. Note that this corresponds to a three dimensional version 
of the expression (25) in [36]. After modifying the package AMBRE to be applicable to three 
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dimensions, a direct application of these packages produces the numerical value of the desired 
function 

7i = /5(2-(0,3-2a),2-(o,3-2(o, 1) = ^ [l+a2e + 0(£2)] (B.2) 

= _2^1^^2Zii:[i + (ii.0365...)e + O(£2)] 

From this, one may easily read tti = — ^ and a2 = 13 + \|/(^). For the integral J2, we have 
found the following numerical value 

72 = 75(1, 3-2(0,1, 3-2(0,2-(o) = pi [l + p2e + 0(£2)] (B.3) 
= -(12.5664...) [l + (7.42278...)£ + 0(£2^ 



This is consistent with the expressions Pi = —An and P2 = 8 + \|/(1). The last integral is evalu- 
ated as 

73 = 75(1,1,1,1,1-0))= [1 + (0.422784... )£ + 6>(£^)] (B.4) 

Alt 

= -[l + (l+X|/(l))£ + 0(£2)]. 
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